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Abstract  - A mixed  element  approach  using 
right  hexahedral  elements  and  right  prism 
elements  for  the  finite  element-boundary 
integral  method  is  presented  and  discussed 
for  the  study  of  planar  cavity-backed 
antennas.  The  mixed  element  method  is 
shown  to  decrease  the  required 
computation  time  for  geometrically 
constrained  geometries  by  reducing  the 
unknown  count  on  the  open  aperture  on  the 
cavity.  By  reducing  the  unknown  count  on 
the  surface,  the  memory  and  computational 
cost  associated  with  the  boundary  integral 
portion  of  the  solution  is  decreased  versus 
solutions  using  only  prism  elements.  The 
accuracy  of  the  mixed  element  approach  is 
shown  to  be  comparable  with  that  of  a 
single  element  approach,  especially  for  far 
field  parameters  such  as  radiation  pattern 
and  radar  cross  section. 


I.  INTRODUCTION 

Efficient  numerical  modeling  of  antennas  is  an 
integral  part  of  the  antenna  design  process. 
Numerical  modeling  can  aid  in  the  rapid 
design  of  an  antenna  prior  to  prototype 
fabrication,  therefore  drastically  reducing  the 
design  time  and  reducing  cost.  One  of  the 
challenges  inherent  in  the  numerical  modeling 
of  antennas  is  how  to  model  the  antenna  in  the 
most  accurate,  yet  efficient  manner. 

Various  numerical  techniques  have  been  used 
to  model  antennas  with  size  on  the  order  of  a 
few  wavelengths.  One  of  the  most  popular  of 


these  techniques  is  formulated  using  an 
integral  equation  (IE)  and  implemented  using 
the  Method  of  Moments  (MoM)  or  one  of  its 
variants.  One  particularly  relevant  example,  a 
spiral  antenna,  is  given  in  [0].  While  this 
method  is  highly  accurate  and  well  studied,  it 
suffers  from  the  fact  that  it  produces  fully 
dense  matrices  leading  to  memory  demands  of 
o(n^)  and  computational  complexity  of 
o(n^)  where  N is  the  order  of  system.  In 
addition,  the  most  efficient  MoM  formulations 
are  based  on  surface  equivalence  and  hence 
are  restricted  to  piecewise  homogeneous 
materials. 

Recently,  techniques  have  been  developed  to 
reduce  the  computational  complexity  of 
integral  equation  formulations  to  ©(Nlogj  N) 
by  exploiting  the  fact  that  many  of  the 
unknowns  are  physically  distant  from  other 
unknowns  [0].  A different  approach,  the 
Finite  Element  (FE)  method  [0-0],  is  based  on 
a partial  differential  equation  (PDE)  approach 
and  therefore  leads  to  very  sparse  system 
matrices  that  can  be  stored  and  solved  in  a 
very  efficient  manner.  In  addition,  since  it  is  a 
PDE-based  approach,  the  FE  method  readily 
permits  analysis  of  inhomogeneous  materials 
in  an  antenna  design.  However,  the  FE 
method  does  not  enforce  the  Sommerfeld 
radiation  condition  for  electromagnetic  waves 
as  part  of  its  formulation  and  hence  is 
susceptible  to  spurious  reflections  from  the 
mesh  truncation  surface.  To  solve  this 
problem,  often  local  conditions,  such  as  an 
absorbing  boundary  condition  (ABC)  or 
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Perfectly  Matched  Layers  (PML)  are  used. 
These  methods  are  particularly  useful  for 
scattering  calculations  since  the  Radar  Cross 
Section  (RCS)  is  a far-zone  quantity  and  hence 
the  effect  of  local  errors  on  the  solution  tends 
to  be  diminished  via  integration.  Solution 
accuracy  is  a particular  concern  for  antenna 
modeling  using  such  an  approximate  condition 
since  important  antenna  parameters  (e.g.  input 
impedance,  mutual  impedance,  etc.)  depend  on 
accurate  local  field  solutions.  A technique 
combining  the  FE  method  and  an  integral 
equation,  termed  the  Finite  Element-Boundary 
Integral  (FE-BI)  method,  is  an  attractive 
alternative  since  it  implements  an  exact 
relationship  between  tangential  electric  and 
magnetic  fields  on  the  mesh  boundary  as  well 
as  the  Sommerfeld  radiation  condition  [0-0] 
via  a properly  constructed  Green’s  function.  A 
third  approach  uses  a harmonic  expansion  of 
the  exterior  field  as  a mesh  closure  condition 
[0].  An  excellent  summary  of  these  various 
conditions  is  given  in  [0]. 

However,  the  FE-BI  method  suffers  from  a 
well-known  drawback:  The  boundary  integral 
portion  of  the  system  takes  up  the  majority  of 
memory  and  time  for  solution.  In  two 
dimensional  cases,  the  use  of  triangular 
elements  on  the  surface  of  the  geometry  to 
achieve  the  greatest  modeling  flexibility  can 
lead  to  electrical  oversampling  which  causes 
larger  than  needed  memory  demand  and 
computational  time.  Quadrilateral  elements, 
while  not  as  flexible  in  modeling,  can  be  used 
to  reduce  the  unknown  count  for  some  cases 
[0],  e.g.  narrow  slots.  The  mixture  of 
triangular  and  quadrilateral  elements  can 
therefore  allow  a flexible  modeling  solution 
while  reducing  the  effects  of  oversampling. 
For  the  case  of  three-dimensional  cavity- 
hacked  antennas,  quadrilateral  elements  can  be 
extruded  into  right  hexahedral  elements  while 
triangular  elements  can  be  extruded  into  right 
prism  elements.  Hence,  the  volumetric 
efficiency  of  quadrilateral  elements  can  be 


combined  with  the  flexibility  of  prisms  to 
reduce  overall  complexity  of  a solution. 

This  paper  presents  the  mixed  element 
formulation  as  applied  to  three-dimensional 
cavity-backed  antennas.  Section  II  of  the 
paper  develops  the  theory  for  the  mixed 
element  formulation.  Section  III  shows  some 
numerical  results  comparing  the  mixed 
element  formulation  with  the  single  element 
formulation.  Section  IV  compares  the 
computational  time  required  for  the  two 
methods.  Section  V presents  some 
conclusions  and  future  directions  for  this 
work. 

n.  FORMULATION 

Figure  1 illustrates  a cavity-backed  aperture 
lying  in  an  infinite,  metallic  groundplane.  The 
computational,  and  antenna,  volume  is  denoted 
by  V while  the  apertiue  is  denoted  by  S.  The 
portion  of  the  groundplane,  including  any 
metal  in  the  aperture,  is  indicated  by  Sp^^ . The 

fields  within  the  computational  volume  are 
denoted  by  (e““  ,H“‘ ) while  the  external  fields 
are  indicated  by  (e“*,H^’“).  The  material 
within  the  cavity  may  in  general  be 
inhomogeneous  and  anisotropic;  however,  for 
this  work  a simplified  formulation  assuming 
isotropic  materials  is  presented.  Expressions 
for  anisotropic  cases  can  be  found  for  prisms 
in  [0]  while  expressions  for  hexahedral 
elements  may  be  similarly  developed. 


Figure  1:  Cavity  backed  aperture  recessed  in  an 
infinite  ground  plane  with  multiple  loading  layers. 
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This  type  of  geometry  represents  flush- 
mounted  cavity-backed  antennas  where  the 
local  surface  in  the  vicinity  of  the  antenna  is 
planar.  Antennas  found  in  this  type  of 
configuration  include  both  narrow  bandwidth 
apertures  (e.g.  microstrip  and  slot  antennas)  as 
well  as  wide  bandwidth  apertures  (e.g.  packed 
elements  such  as  I-dipoles  as  well  as  spiral  and 
log-periodic  antermas).  A systematic  modeling 
of  these  antennas  through  the  use  of  highly 
efficient  computational  algorithms  is  important 
to  determine  the  optimal  configuration  of  these 
systems. 

The  problem  domain  is  separated  into  two 
distinct  regions,  the  region  composed  of  the 
interior  volume  of  the  geometry,  and  the 
region  composed  of  the  surrounding  free  space 
area  above  the  ground  plane.  The  interior,  FE, 
region  is  then  coupled  to  the  exterior,  BI, 
region  by  enforcing  tangential  magnetic  field 
continuity,  nxH®**  =hxH“‘,  on  the  aperture 
surface  (S')  while  tangential  electric  field 

continuity,  h x E“‘  = n x E“‘ , on  that  same 
surface  is  assured  via  the  use  of  similar  basis 
function  for  the  interior  and  exterior 
representations.  The  final  hybrid  expression, 
usually  known  as  the  finite  element-boundary 
integral  equation,  is  given  by 

-!-£(vxi:,).(vxE“)  dr- 

dr  (1) 

where  Lj  is  the  vector  testing  function 

associated  with  the  i***  row  while  J”*’’ 
represents  an  impressed  current  source  as  the 
antenna  feed.  The  relative  material  parameters, 
e^and  p,,  are  associated  with  each  element 
and  may  vary  on  an  element-by-element  basis 
while  kg  and  are  the  firee-space 

wavenumber  and  wave  impedance, 
respectively.  The  electric  dyadic  Green’s 


function  of  the  second  kind,  > is  taken  to 

be  the  half-space  Green’s  function  [0]. 

The  unknown  electric  field  is  expanded 
using  the  same  functions  as  used  for  testing 
(e.g.  Galerkin’s  method)  is 

(2) 

j=l 


where  Ej  indicates  complex-valued  unknown 


expansion  coefficients.  Inserting  (2)  into  (1), 
the  final  expression  is  given  as 


^((vxLHvxl)  dr-«j£, |;i,  ^ dr 

O'  ' ' 

£_^-(zx4(:«;yK,y)xz)-L,.  dSW 


(3) 


This  matrix  equation  is  separated  into  two 
parts,  one  representing  the  finite  element 
portion,  written  as  I™,  and  the  other 

representing  the  boundary  integral  portion, 
written  as  l“  as 


=-yL(^xL,)  (VxL,)  dV 

dv  , (4) 

I?= 

IJ 

“ Is,  4 ^ ^e2 (x,  y I x',  y') X z)-  Lj  dS'dS 


(5) 

where  the  latter  has  support  only  when  both 
test  and  source  edges  lie  in  the  aperture.  These 
integrals  represent  the  matrix  entries  in  the 
following  linear  system 

FE  , tBI 


I 


y 

tFE 

^ij 


-fl- 
y 


tFE 


I 


FE 


(6) 


where  E®  indicates  unknown  field  expansion 
coefficients  associates  with  the  aperture 
surface,  E“*  are  the  expansion  coefficients  for 

the  interior  basis  vectors  and  f is  the 
excitation  term  associated  with  an  interior 
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current  source.  Evaluation  of  (4)  and  (5) 
require  specification  of  the  basis  functions, 
Lj.  These  will  be  different  for  the  different 

types  of  elements  used  in  the  mesh;  hence  their 
will  be  two  different  specifications  of  ly® , one 
for  prisms  and  one  for  hexahedra,  and  four 
different  specifications  of  I?  to  represent  all 

possible  combinations  of  surface  element 
coupling.  The  functions  for  a right  prism  are 
presented  first  followed  by  those  for  a right 
hexahedral. 

Right  Prisms 

Right  prism  elements  have  been  used 
successfully  in  the  FE-BI  method  in  the  past 
[0].  These  basis  functions  posses  the  required 
properties  to  be  used  in  the  FE-BI  method. 
They  enforce  tangential  field  continuity  across 
element  faces  and  are  curl  conforming.  The 
elements  used  in  this  work  are  also 

divergence-free  and  CT/LN  elements  [0]. 

Prisms  have  the  advantage  of  providing  great 
flexibility  in  modeling  geometries  that  are 
irregular  in  two  dimensions  but  regular  in  the 
third  dimension,  such  as  cavities  recessed  in 
ground  planes,  the  case  presented  here.  Prisms 
also  have  the  advantage  of  making  it  quite 
simple  to  extrude  a three-dimensional  mesh 
from  a two-dimensional  mesh  composed  of 
triangles.  In  addition,  the  prism  basis 

functions  are  derived  from  the  two 

dimensional  basis  functions  for  triangles  so 
they  reduce  to  the  two-dimensional  triangular 
basis  functions  used  in  the  boundary  integral; 
hence,  prisms  are  similar  to  triangle  elements 
and  therefore  automatically  enforce  the 

essential  boundary  condition  in  the  aperture. 
However,  for  some  geometries  such  as  narrow 
slot  spiral  antennas,  prism  elements  have  the 


disadvantage  of  oversampling  the  aperture  and 
therefore  inefficiently  using  resources.  Figure 
2 illustrates  one  possible  realization  of  a right 
prism  element  where  the  local  node  numbering 
scheme  is  shown  along  with  the  encircled  local 
edge  numbering  scheme. 


Figure  2:  Right  prism  element  shown  with  its  defined 
local  nodes  and  local  edges  (encircled  numbers). 

The  nodes  are  ordered  in  a counter-clockwise 
direction  in  order  to  ensure  that  the  normal 
vector  of  the  element  points  towards  the  top  of 
the  element  and  therefore  points  out  of  the 
computational  volume  in  the  aperture  of  a 
cavity-backed  antenna.  The  edge-based 
expansion  functions  for  the  prism  element  are 
derived  from  the  Rao-Wilton-Glisson  (RWG) 
basis  functions  for  triangles  [0]  with  a linear 
depth  variation  to  allow  for  three-dimensional 
support.  Assigning  prism  vector  expansion 
functions  asL^  = , the  nine  vector  functions 

associated  with  the  prism  are  given  by 
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W;,,-  = Y =^^^||r[(x-x,)y-(:>;-);,)i]  , 

W^,  =M,  =^^^^[(x-x,)y-(y-y,)i]  , 

'^zi  [(^H  ) + {yn  - yk2 ) ^ + (^i-2  - ) y ] , 


Z = h2,3, 
Z = 4,5,6, 
X = 1A9, 


where  /;is  the  length  of  the  i"'  edge,  is  the 

sign  of  the  edge  to  ensure  field  continuity 
between  elements  , (x;,y^)  are  the  global 

coordinates  the  local  nodes,  (z/,z„)  are  the 
global  coordinates  of  the  upper  and  lower 
faces  of  the  prism,  c = - z,  is  the  height  of 

the  prism,  and  is  the  area  of  the  triangle 
that  forms  the  top  (or  bottom)  of  the  prism.  For 
the  vertical  edges,  and  k2  are  defined  in 
Table  1 

Table  1.  Definition  of  the  indices  used  to  construct 
the  vertical  prism  expansion  functions. 


JL 

ki 

k2 

7 

2 

3 

8 

3 

1 

9 

1 

2 

where  % is  the  local  edge  number  and  (kj,lc2) 
indicate  local  nodes.  The  curls  of  these  basis 
functions  are  defined  as 

V X V;.  = - x.)\  + {y-  y.)y  - 2(z  - z,  )z], 

^ ^ + (y~  + 2(z„  - z)z], 

2cS 

[fe-%)x+(yH-y^2)y] 

'■  2S^ 

(8) 

Note  that  these  expansion  functions  are 
functionally  identical  to  [0];  however,  since 
they  are  expressed  in  global  rather  than  local 
coordinates,  it  is  relatively  easy  to  use  these 
functions  in  conjunction  with  anisotropic 
materials  specified  in  terms  of  global 
properties. 


(7) 

Right  Distorted  Hexahedrals 
A common  example  of  a regular  hexahedra 
element  is  the  brick  element  where  all  the 
edges  are  either  parallel  or  orthogonal  to  any 
other  edge  in  the  element  [0].  Distorted 
hexahedral  elements  have  edges  that  are  not 
necessarily  parallel  or  orthogonal  to  the  other 
edges.  Distorted  hexahedral  have  been  used  by 
[0]  and  discussed  in  [0];  however,  not  in 
conjunction  with  other  elements.  The  edges 
in  the  extrusion  direction  are  orthogonal  to  the 
quadrilateral  element  used  to  form  the 
hexahedron.  Distorted  hexahedral  elements  are 
important  because  they  have  the  ability  to 
model  irregular  surface  geometries  with 
potentially  fewer  edges  than  prisms  (e.g. 
narrow  slots).  This  leads  to  less  computational 
and  memory  demand  for  the  same  geometry 
modeled  with  hexahedral  elements  as  opposed 
to  prism  elements.  Distorted  hexahedral 
elements,  however,  have  the  disadvantage  of 
not  having  closed-form  matrix  entry  formulae 
and  hence  require  more  computational  effort 
as  compared  to  prisms  with  similar  field 
representation  capability. 

The  basis  functions  used  for  hexahedral 
elements  are  the  so-called  rooftop  functions. 
Rooftop  basis  functions  are  most  often  used  in 
brick  elements  [0-0].  Brick  elements  are  very 
easy  to  use  but  suffer  from  the  fact  that  they 
can  only  model  Cartesian-t5q)e  geometries 
effectively.  In  order  to  model  more  irregular 
geometries  distorted  hexahedral  elements  can 
be  used.  Due  to  their  distortion  however  the 
integrals  in  (5)  can  be  very  difficult  to 
compute  numerically.  Numerical  integration 
over  brick  shaped  volumes  is  readily 
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implemented  and  for  this  reason  it  is  mapped  into  a brick  before  integration, 
advantageous  if  a hexahedral  element  is 

Figure  3 shows  the  transformation  of  a distorted  hexahedral  element  in  the  {x,y,z)  coordinate 
system  to  a unit  cube  in  the  coordinate  system. 


Figure  3:  Hexahedral  in  (x,  y,  z)  coordinates  mapped  into  a cube  in  (^,  i|,  Q coordinates. 


where  the  edges  of  the  hexahedral  and  the 
brick  are  defined  in  Table  2. 

Table  2 : Local  edge  numbering  for  a distorted 
hexahedral  element. 


Edge 

Nodel 

Node2 

1 

1 

2 

2 

3 

4 

3 

1 

4 

4 

2 

3 

5 

5 

6 

6 

7 

8 

7 

5 . 1 

8 

8 

6 

7 

9 

1 

5 

10 

2 

6 

11 

3 

7 

12 

2 

6 

From  [0]  the  vector  edge-based  basis  fimctions 
can  be  written  as 

n,  = %(1  + /;,7)(i  + ^,4')V| 


for  the  edges  parallel  to  the  ^ direction 

N,  = ^(l+^,.^)(l  + C.4')Vt7 

(10) 

for  the  edges  parallel  to  the  p direction,  and 
finally 

N,  = ^(1  + ^,^)(1  + /7,.;7)V4' 

O 

(11) 

for  the  edges  parallel  to  the  ^ direction,  where 
/,  denotes  the  length  of  the  i"'  edge.  Based  on 

the  definition  of  the  gradient,  these  equations 
can  be  equivalently  written  as 

N, 

N,=i^{i+^,^Xi+C,C)n, 
Ni=^(l  + n,riXl  + a)C- 


8 


(9) 


(12) 
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Therefore,  the  basis  functions  for  hexahedral 
elements  are  now  defined  in  the  mapped 
coordinate  system  where  numerical  integration 
can  be  performed  over  a cube  rather  than  a 
hexahedral  volume. 

The  integrands  in  (5)  are  still  represented  in 
terms  of  global  Cartesian  coordinates. 
Defining  the  Jacobian  matrix  as 


dx 

dy 

dz 

d4 

d^ 

dx 

dy 

dz 

d 1] 

d T] 

d T] 

dx 

dy 

dz 

d^ 

\ 

(13) 

the  elemental  volume  of  the  integral  can  be 
written  as  [0] 


dV  = dxdydz  = det[J]d^d?jd^  . 

(14) 

From  [0]  the  following  vector  operations  can 
be  defined 


z-(VxN)  = 


det[J] 


dN”  dN^ 
dr] 


(VxlS)=- 


det[J] 

+n 

\dn  dcj  \ 


dc 


(15) 

/av"  av^^l 


drj 


The  above  equations  can  be  used  to  derive  the 
equations  needed  to  implement  the  FE-BI 
method  using  hexahedral  elements. 

Mixed  Element  Formulation 
Traditional  formulations  of  the  FE-BI  method 
use  one  of  the  previous  two  elements 
described  as  the  sole  type  of  element  used  to 
model  the  geometry.  This  works  well  in 
practice  and  many  implementations  of  this 
type  have  been  successful.  Each  method  has 
its  advantages  and  disadvantages  as  described 
previously.  To  gain  benefits  not  obtainable 


using  one  element  or  the  other,  the  two  types 
of  elements  can  be  combined.  This  allows  the 
use  of  prism  elements  where  their  flexibility  is 
needed,  such  as  at  sharp  contours  and  areas  of 
rapid  varying  fields,  while  still  allowing  the 
use  of  hexahedral  elements  where  fewer 
unknowns  are  needed,  such  as  areas  of  slowly 
varying  fields  or  for  narrow  slots.  Only  a few 
additional  equations  need  to  be  introduced  in 
order  to  make  this  mixed  element  formulation 
possible.  To  be  more  specific,  only  the 
boundary  integral  terms  where  a prism  edge  is 
interacting  with  a hexahedral  edge  need  to  be 
derived. 


The  equations  for  mixed  element  interactions 
via  the  boundary  integral  can  be  derived  by 
using  the  basis  functions  for  the  triangular  and 
quadrilateral  elements  in  (5).  For  the  case 
when  triangles  are  used  as  test  elements  and 
quadrilaterals  are  used  as  source  elements,  the 
following  equations  are  used 

jm  I j J[Mf-(^+^)-N,]^det[J,]^^3&7dS', 

^ s-\-\  ^ 

(17) 

(18) 

where  R = -^{x  - x')^  + (y  - y')^ 

(16) 


Conversely  when  quadrilaterals  are  used  as 
test  elements  and  triangles  are  used  as  source 
elements,  the  following  equations  are  derived. 

Ip-  ^ ^ 

jm  J j Jps;.  ■(Sx-F^).\^.]^dd[J,]riS'^/4i/7, 


1 vv 


(19) 


,-J^ 


(20) 


Using  the  following  relations 


V-(zxN,)  = — 

det[J,] 


dNf  dNf 

d7] 


,(21) 
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V -(Z  X W ,)  = - 

the  equations  needed  to  implement  the  mixed  element  formulation  are  written 


(22) 
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where  p denotes  the  fact  the  either  the  test  or 
source  edge  are  associated  with  a prism 
element. 


Figure  4:  Radar  Cross  Section  of  a 6cm  x 5cm  x 2cm 
cavity  with  a 3cm  x 2cm  slot  aperture,  00-pol. 


EQ.  Numerical  Results 

To  validate  this  mixed  element  formulation, 
two  simple  test  cases  were  constructed.  The 
first  was  a 3cm  x 2cm  slot  antenna  cut  into  a 
6cm  X 5cm  X 2cm  cavity.  The  cavity  was 
filled  with  a dielectric  material  with  a 
dielectric  constant  of  = 2.17 . The  geometry 
was  modeled  once  with  only  prism  elements 
and  a second  time  with  a mixture  of  prism 
elements  and  hexahedral  elements.  The  RCS 
of  this  antenna  is  shown  in  Figures  4 and  5. 


Figure  5:  Radar  Cross  Section  of  a 6cm  x 5cm  x 2cm 
cavity  with  a 3cm  x 2cm  slot  aperture,  (p^-pol. 
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The  results  above  show  that  the  mixed  element 
formulation  is  nearly  identical  to  the  single 
element  formulation  when  modeling  the  radar 
cross-section  of  the  slot  antenna. 

The  second  test  case  consisted  of  a 3cm  x 2cm 
patch  antenna  residing  in  a 6cm  x 4cm  x 
0.0762cm  cavity  filled  with  a dielectric 
material  with  Sr  = 3.2.  The  patch  antenna  was 
excited  by  a probe  feed  located  at  x = 3cm,  y = 
2.1cm.  As  with  the  slot  antenna  this  antenna 
was  modeled  once  with  only  prisms  and  again 
with  a mixture  of  prism  elements  and 
hexahedral  elements.  Figure  6 shows  the 
radiation  pattern  of  the  patch  antenna  at  5.4 
GHz  and  Figure  7 shows  the  input  resistance 
of  the  probe  feed  from  5-6  GHz. 


Figure  6:  Normalized  radiation  pattern  for  a 3cm  x 
2cm  patch  antenna  in  6cm  x 4cm  x 0.0762cm  cavity. 


Figure  7:  Input  resistance  of  a 3cm  x 2cm  patch 
antenna  in  6cm  x 4cm  x 0.0762cm  cavity. 


The  above  results  show  that  the  mixed  element 
formulation  again  matches  the  prism-only 
formulation  results  very  closely.  Different 
mesh  densities  around  the  feed  point  for  the 
two  different  methods  cause  the  discrepancy  in 
the  value  of  the  input  resistance  at  the  resonant 
point  of  the  patch  antenna.  Due  to  the  extreme 
local  nature  of  input  impedance  small  changes 
in  electric  field  values  caused  by  slightly 
different  sampling  rates  at  the  feed  point  can 
lead  to  moderately  different  results  for  the 
value  of  the  input  impedance.  Also,  probe 
feeds  such  as  the  one  used  in  the  above 
example  can  be  extremely  sensitive  to  small 
changes  in  the  field,  which  makes  them 
somewhat  unreliable  for  accurate  simulations 
of  input  impedance  value.  The  most  important 
result  in  the  above  example  is  that  the  structure 
has  a resonant  point  at  the  same  frequency 
using  both  the  prism  element  formulation  and 
the  mixed  element  formulation. 

rv.  Solution  Efficiency  Comparison 

To  compare  and  contrast  the  computational 
demand  required  to  model  geometries  with 
only  prism  elements  as  opposed  to  a mixture 
of  prism  and  hexahedral  elements  a few,  more 
complex,  cases  were  considered.  These  cases 
were  chosen  since  they  represent  a class  of 
problems  where  the  triangles  used  with  prisms 
to  represent  the  aperture,  “oversample”  the 
aperture  from  an  electromagnetic  viewpoint. 
This  is  best  understood  by  considering  a 
narrow  slot  antenna.  Since  the  eleetrie  fields  in 
the  slot  have  only  a component  perpendicular 
to  the  slot  sides,  quadrilateral  elements  very 
efficiently  model  this  slot.  In  contrast,  using  an 
identical  sampling  density,  triangle  elements 
would  require  one  additional  degree-of- 
freedom  (e.g.  edge)  per  equivalent 
quadrilateral  element. 

The  first  of  these  cases  included  a four-arm 
spiral  antenna  with  one,  two,  and  three  turns. 
The  second  was  an  I-dipole  array  consisting  of 
one,  four,  nine,  and  sixteen  dipoles.  The  time 
required  computing  a single  radiation  pattern 
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cut,  with  an  INTEL  XEON  550  MHz 
processor  running  the  Linux  operating  system, 
was  recorded  for  each  of  the  aforementioned 
geometries.  Figure  8 shows  examples  of  the 
different  geometries  used  for  the  results  along 
with  the  meshing  used.  Table  3 describes  the 
different  mesh  parameters  for  the  example 
geometries.  A single  pattern  cut  was 
computed  in  the  X-Z  plane  at  a frequency  of 


12  GHz  for  the  spiral  antenna  and  4GHz  for 
the  I-dipole  array.  Figure  9 illustrates  the  time 
required  to  compute  a radiation  pattern  cut  for 
the  four-arm  spiral  antenna  versus  the  number 
of  turns  in  the  spiral.  Figure  10  illustrates  the 
time  required  to  compute  a single  radiation 
pattern  cut  for  the  I-dipole  array  versus  the 
number  of  dipoles  in  the  array. 


(c)  (d) 

Figure  8:  A four-arm,  one-tum  spiral  and  a 3x3 1-Dipole  array  (a)  spiral  with  mixed  elements,  (b)  spiral  with 
prism  elements,  (c)  I-Dipole  array  with  mixed  elements,  (d)  I-Dipole  with  prism  elements. 


Table  3:  Mesh  comparison  for  a four-arm,  one-turn  spiral  antenna  and  a 3x3 1-Dipole  array. 


Spiral  Mixed 

Spiral  Prism 

I-Dipole  Mixed 

I-Dipole  Prism 

Surface  Elements 

1313 

2338 

2280 

3600 

Surface  Edges 

2544 

3569 

3980 

5437 

Surface  Unknowns 

1296 

1941 

2880 

4013 

Total  Elements 

14443 

25718 

15960 

25200 

Total  Edges 

44080 

56380 

43747 

56362 

Total  Unknowns 

37684 

48579 

37705 

48539 
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Figure  9:  Required  computation  time  for  a single 
radiation  pattern  cut  for  a four  arm  spiral  antenna. 

X 10“ 


F^ure  10:  Required  computation  time  for  a single 
radiation  pattern  cut  for  an  l-dipole  array  antenna. 


As  can  be  seen  from  the  above  results  the 
mixed  element  formulation  scales  much  better 
for  geometrically  constrained  problems,  such 
as  spiral  antennas  and  I-dipole  arrays,  as 
compared  to  the  prism-only  formulation. 

V.  Conclusions 

A mixed  element  formulation  was 
developed  to  model  three-dimensional  planar 
cavity-backed  aperture  antennas.  This 
formulation  was  shown  to  be  equivalent  to  a 
formulation  using  only  prism  elements  but  was 
more  efficient  in  terms  of  memory  and 


computational  demand.  The  main  memory 
and  computational  savings  were  due  to  the  fact 
that  the  mixed  element  formulation  produced 
fewer  surface  unknowns,  for  a comparable 
electromagnetic  representation,  than  the  prism 
element  formulation  for  the  same  geometry. 

Slight  discrepancies  between  the  mixed 
element  and  prism  element  formulations  were 
seen  in  the  modeling  of  impedance  for 
complex  geometries.  Dissimilar  mesh  densities 
surrounding  the  feed  points  causing  different 
input  powers  to  be  produced  at  the  feeds 
caused  these  discrepancies.  These 
discrepancies  were  not  seen  in  radar  cross- 
section  simulations  as  RCS  measurements  do 
not  include  a feed  model. 
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